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ABSTRACT 


Pcase (1987) related the effects of atmospheric forcing. mainly temperature and wind 
fields, to the size of coastal polvnas. Using Pease's formulation and Killworth's (1977) 
plume model as applied bv Melling and Lewis (1982), the effects of atmospheric forcing 
on brine injection into the Arctic pvenoeline are investigated. This paper will discuss the 
likelihood of coastal polvnas as a source for denser abyssal waters. 

A standard case was developed for the model with initial conditions taken from 
Melling and Lewis (1982) and Pease (1987) for comparison with individual sensitivity 
experiments. Ten environmental parameters were individually examined for their influ- 
ence on the plume depth after 90 days. The standard case resulted in a 90-day plume 
depth of 436 meters. A submarine canyon case was simulated, resulting in plume pene- 
tration to over 1300 meters in 90 davs. Further experiments used actual T-S soundings 
from Aagaard et al. (1981) and Ostlund et al. (1987). Finally, a 20 kilometer wide plume 


is shown to penetrate to almost 600 meters in 90 days. 
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I. INTRODUCTION 


ITI TTT remotelv located part ol the world's ocean. has de- 
manded the attention of those who wish to understand the Earth's climate. The renewed 
emphasis in recent vears given to understanding the Earth's climate has underscored the 
RucEolehnevArcte Occanum tbe scheme of the global'energv balance. The Arctic Ocean 
is also a key player in the world oceans’ thermohaline circulation. Because of the im- 
portance of the Arctic Ocean in the global picture, understanding the nuances of the 
Arctic Ocean's physical oceanography and its sensitivities to environmental change is 
essential before climatologists, oceanographers, and meteorologists can develop a more 
realistic and longer range global climatological concept. Two puzzling aspects of the 
Arctic Ocean are its isothermal halocline and the more saline deep waters. How is the 
halocline maintained in a nearly horizontally independent state? What is the mechanism 
behind the formation of the Arctic Ocean's deep waters? These topics will be discussed 
in this thesis. 

Throughout the Arctic Ocean, the characteristics of the upper lavers are remarkably 
independent of the horizontal location of the sample. The water column consists of a 
mixed laver, halocline, thermocline, and deeper waters. The mixed laver 1s a relatively 
fresh water layer approximately 30 meters deep. River run-off and the summer cvcle of 
melting 1ce are considered the main contributors to the fresh water laver. The halocline 
ne laver ol water beneath the mixed laver. In this layer, salinities rise from 31 to as 
high as 34.6 ppt at the 200 meter depth level (Killworth and Smith, 1984). The halocline 
IS, IN essence, an isothermal laver with temperatures just above the freezing point. The 
thermocline 1s a layer below the halocline, in which temperatures rise to approximately 
0* C at 250 meters (Melling and Lewis,1982) and continues to rise to 0.5* C at 400 meters 
due to the laver of warm Atlantic water. Below this depth, temperatures decrease with 
increasing depth. The deep waters of the Arctic Ocean are among the most dense waters 
found in the world's oceans and have been suggested as being an additional source for 
Ventilating the world’s oceans’ deep thermohaline circulation (Aagaard,1981). Although 
the physical oceanography of the Arctic Ocean is well described in literature, the exact 
mechanism(s) for the maintenance of the halocline and deep water formation are not 


fully understood. 


The halocline is generally isothermal except in the Canadian Basin, where warin 
Pacific Waters entering via the Bering Strait causes a warm laver at about the 75-100 
meter depths, and presents a paradox when examining its origin. Sandwiched between 
a cold. relatively fresh laver and a warm. saline laver. an obvious solution to the main- 
tenance of the halocline might be simple vertical mixing between the two lao o 
ever, vertical mixing of the cold, less saline surface water with the warmer. more saline 
Auantic water 1s not considered to be a source since heat is not mixed along with salt. 
Aagaard et al. (1981) showed a discontinuity existed in the T-S correlation slope at 
about 150 meters. Either salimizationscr surface waters by freezing (1.e. brine rejeerens 
or cooling of Atlantic water could be two possible mechanisms to produce halocline 
water. The authors (Aagaard. et al., 1981) demonstrate mixing between the two layers 
cannot be a source of halocline water by examining the required salinity of the surface 
water. Using a T-S diagram, and assuming linear mixing at a rauo of 3:2 (surface water 
to Atlantic water), the surface water would need to have a salinity of 34.7 ppt. Salinity 
this high has not been observed in the Arctic surface laver. Therefore, the maintenance 
of the Arctic halocline 1s generally thought to be attributable to lateral advection of cold, 
saline water from the continental shelves. The shelf waters are made more saline due to 
brine rejection associated with ice growth. Aagaard (1981) proposed that continental 
shelf areas where an annual growth of ice less than 2 meters is necessary to modifv 
summertime salinity profiles to wintertime profiles may be the source areas of the 
halocline waters. Upwelling of Atlantic water onto the shelves and subsequent cooling 
has also been suggested as a possible source. Upwelling along submarine canyons has 
been reported as early as 1929 (Sverdrup,1929). As further evidence eluding quomm 
halocline's possible shelfwater origin, water on the continental shelves has been shown 
to have temperature and salinity characteristics similar to the waters found offshore at 
a depth of 150 meters (Melling and Lewis, 1982). 

In their paper, Melling and Lewis (1982) discuss the lateral advection of low tem- 
perature high salinity shelfwaters into the halocline, with particular attention paid to the 
effects of brine rejection due to the growth of ice. Aagaard (1981) used two methods 
of computing the necessary production rate of high salinity shelfwater to maintain the 
halocline. By considering the production rate required to renew the halocline in 10 vears 
and necessary mixing ratios to achieve desired T-S characteristics, Aagaard suggests a 
2.5x10% to 4.0.x10%m — sec"! for a production rate of cold saline shelfwater, assuming a 


halocline renewal period of 10 vears. In order to increase the ambient shelfwater salinity 


tv 


to values which would form sufficiently dense water to descend the slope into the 
halocline, Melling and Lewis (1982) emphasize the important role of ice dynamics. spe- 
cifically the role of divergence. Areas of large wintertime ice cover divergence allow 
rapid ice growth hence increased salinization of the underlying water column. Leads and 
polvnas are locations of rapid ice growth due to the divergence of the insulating ice 
cover. 

cnn acerco s7s5) examined the degree of salinization of the Underlying water 
column beneath a polvna. In their study. measurements of salinity and water current 
velocity beneath a polvna located at St. Lawrence Island indicated polvnas to be an area 
ENNIebuceerowthorates resultine im densitv driven currents due to increased salinization 
E ehtimmacher eral, 1955). Coupling the salt source. 1.¢. polyna, with the interleaving 
plumes of dense shelfwater found by Melling and Lewis (1982) may provide a possible 
mechanism for the maintenance of the Arctic halocline. The polvna provides an atmo- 
spheric dependent source of salt into the water column possiblv forming plumes which 
descend down into the halocline, as described bv Melling and Lewis (1982). Along with 
the mechanism associated with the haloclines maintenance, the formation of the 
Canadian Basin's deep water has also been attributed to having a source of low tem- 
perature, high salinity shelfwater. 

The dense waters which ventilate the world’s oceans flow out of the Arctic Ocean 
over the Greenland-Scotland ridge system at a depth of about 2600 meters. A single 
source was advocated by Helland-Hansen and Nansen (1909) in which deep water is 
formed in the Greenland Sea. Nansen (1906) had written earlier of the possibility of the 
Greenland Sea Deep Water being modified by cooling in the Barents Sea. This theory 
onutted the possibility of deep water formation within the two basins of the Arctic 
Ocean, namely the Eurasian and Canadian Basins. In 1981, Aagaard postulated that 
deep water is also formed within the Arctic Ocean due to the salinity distribution. 
Greenland Sea Deep Water has a salinity of approximately 34.90 ppt, which is the least 
saline of Arctic deep waters. Therefore, at least two sources of deep water formation 
are thought to exist. One source is the traditional Greenland Sea formation area and the 
new source proposed bv Aagaard, 1s within the Arctic Ocean's two basins. 

The characteristics of water mass properties found in the Canadian and Eurasian 
Basins are similar until about the depth of the Lomonosov Ridge. At about 1500 meters 
depth, the Canadian Basin becomes more saline and is slightly warmer than the Eurasian 


Basin Deep Water (Aagaard, 1981). Under the traditional viewpoint of Greenland Sea 


Deep Water filling the two basins through lateral advection, salinities should not exceed 
34.9 ppi. Tlowever, during the LOREX (Lomonosov Ridge l:xperiment; Weber, 1979) 
salinities in the t sins were detcrmincd to be greater than 34.9 ppt. This suggests an 
additional source of salt 1f not a completely new source of deep water. Current meters, 
placed ovcr the Lomonosov Ridge, gave evidence to the existence of some tidal and 
other episodic surges carrying salinized Greenland Sea Deep Water over the ridge into 
the Canadian Basin. Even with the surges over the Lomonosov Ridge into the Canadian 
Basin, water from the Eurasian Basin must have an additional salt source to obtain the 
charactcristics of water found in the Canadian Basin (Aagaard, 1981). Residence times 
can be calculated by comparing the mixing ratios of shelf water to Atlantic Intermediate 
water necessary to achieve the T-S characteristics found in the Canadian Basin Deep 
Water. Considering a small flux of shelfwater (.003 to .06 Sv) necessary to ventilate the 
Canadian Basin over several hundred years may also explain the warmer temperatures 
found in the Canadian Basin (Aagaard ct al., 1985). By using tracers such as ó'*O and 
“C, Ostlund et al. (1987) agree with the magnitudes of the renewel time scales for both 
basins presented by Aagaard et al.(1985). Ostlund et al. (1987) and Aagaard et al. (1985) 
place the residence time of the Eurasian Basin in the order of decades vice several hun- 
dred years (700 years, Ostlund, et al., 1987) for the Canadian Basin. The requirement 
remains of an additional salt source at small fluxes to ventilate the Canadian Basin over 
several hundred years point to an episodic mechanism. Killworth and Smith (1984) write 
that shelfivater sinking at the basin edges may occur but a source of waters sufficiently 
dense to penetrate the Arctic Ocean's stable mid-layer remains a problem. Again, using 
a coastal polvna as a salt source, could shelfwater be modified under extreme conditions 
to produce such dense water which would then descend to the deeper depths via gravity 
flow plumes? 

This tnesis will combine Pease’s numerical formulations relating atmospheric forc- 
ing, mainly wind and air temperature, to the opening and closing of coastal polynas, with 
a streamtube gravity flow model developed by Smith (1975). With simplified 
hydrography and suitable environmental conditions, this model can produce events of 


deep plume penetration along with à more consistent source of high salinity, low 


temperature water for the maintenance of the Arctic Ocean's halocline. Chapter 2 will 


describe the numerical model used in this studv along with an examination into 
drawbacks associated with Smith's streamtube model as described bv other studies. 1.e. 
Bee lee 197. and Melline and Lewis (1982). Chapter » will describe 
results obtained by altering environmental inputs into the model and Chapter 4 will 


discuss the results and implications. 


In 


Il. NUMERICAL METHODS AND PARAMETERS 


A. THE NUMERICAL MODEL 

With the increased desire to understand the nature of the Arctic halocline and its 
origin, there existed sufficient pieces of the puzzle to develop a coupled model. The 3-D 
turbulent entraining plume model developed bv Smith (1975) has been used bv several 
studies ( Killworth (1977), Killworth and Carmack (1978). and Melling and Lewis (1982)) 
to gain insight into gravity flow currents arising from brine rejection into the water col- 
unin due to ice growth. Melling and Lewis (1982), compute the initial densitv of the 
plume by estimating the average ice grown annually in a shallow shelf and and the re- 
sultant salinity contribution. In this study, Pease’s (1987) formulation relating a 
polvna's size to offshore wind speed and air temperature, provides a specific source of 
brine Which is then coupled with Smith's plume model, as used by Melling and Lewis 
(1982), in order to study the possible contribution of the polyna produced denser waters 
to the Arctic halocline and possible formation of Canadian Basin Deep Water. 

1. The coastal polyna 

Pease (1987) describes the polyna’s maximum width as an equilibrium between 

frazil ice production and its advection rate. The frazil ice production rate 1s a function 
of the heat exchange at the air-sea interface. Therefore, with wind speed and al 
perature as independent variables, one can compute the size of the polyna and its frazil 
Ice production rate. 


The polvna size can be found bv solving: 





d 1 
dt ' II, 
where ak Is the change in polvna width with ume, V, is the advection rate of frazil ice, 
( 


X, 1s the polvna width, £F, is the frazil ice production rate, and Z is the collection depth 
of grease ice or nilas. The solution of equation (1) assumes the meteorological condi- 
tions are steady throughout the polvna's opening time, thus allowing treatment as a 


linear differential equation. Furthermore, if the meteorological event 1s assumed to be 


of sufficient duration so that the polyna reaches maximum width. equation (1) can be 


written as: 


mM UN 
Ap naxanum) = ——, (2) 
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This allows solving equation (2) for the maximum width after computing the 
frazil ice production rate, F. In this studv, meteorological events of interest are 4-10 
days duration, therefore this assumption allowing equation (2) is valid. Ou (1988) de- 
scribes a more detailed coastal polyna model and concludes the ice edge is less affected 
bv higher frequency atmospheric disturbances than longer period, svnoptic tvpe vari- 
ations. Although Ou (1988) includes additional phvsics to model the temporal changes 
of the polvna's ice edge, the steady state solution is still that as given by Pease (1987) 
and shown above in equation (2). 


The frazil ice production rate is given by Pease (1987) as: 


[ce, T; S EE p 2 
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where the evaporative heat flux has been neglected due to its small contribution relative 
to the uncertainty of the sensible heat flux. Q,, is the upward longwave radiation, ce,T? 
is the downward longwave radiation, L is the latent heat of freezing for salt water. p, is 
the density of young sea 1ce, and. o,C;, C, V (T, — T.) is the sensible heat flux. 
2. Salinity distribution within the water column 
The salinity contribution to the water column from frazil ice growth is assumed 


to be uniform with depth and is computed using: 


— Lä sz Abr, 
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from Killworth and Smith (1984). S... is the new salinity value of the water column after 
brine rejection. S is the average salinity of the shelfwater, S, is the salinity of the ice, 4, 
is the effective area of salinization, and f qdz is the areal flux of the underlying water. 
Using F, from the polyna model. one can compute the salinity of shelfwater exiting a 
polvna. 
3. Equation of state 
The density of the current 1s computed using the International one atmosphere 


equation of state developed by Millero and Poisson (1981). The equation is given below. 


> 


p=P.+AS+BS=CS (5) 
where, 


A = 8.24493 x 107! — 4.0899 x 1077 + 7.6438 x 101? — 8.2467 x 10 7 + 5.337 MN 
B = —5.72466 x 107? + 1.0227 x 1077 — 1.6546 x 10 — 6r? 


ea, aloe 


and p, 1s the standard densitv of seawater. This allows the calculation of the current's 


relative buovancv using: 


A al (9) 


[n this fashion, the results from the polyna model can be used as inputs to the 
streamtube model. 
4. The streamtube model 

Smith's (1975) model of boundary currents includes the effects of friction, 
entrainment, and Coriolis forces. The 3-D plume was modified by Killworth (1977) by 
altering the entrainment rate. Kıllworth (1977) related the entrainment Tate OS 
product of the plume’s area and velocity vice only the velocity (Smith, 1975). Melling 
and Lewis (1982) used Killworth’s version of the plume model with an entrainment co- 


efficient dependent on the Richardson number. The model used in this paper is the 


Melling and Lewis (1982) adaptation with the inclusion of several stratified lavers to al- 
low a single model initialization. 

Smith’s plume model uses two coordinate systems. The Cartesian system 1s 
aligned with the bottom so the X-axis lies along the shore and the Y-axis lies along the 
slope. The second coordinate system is the curvilinear svstem with 4 and y axes. The 
position of the plume axis can be defined by the value of £. The angle which the plume 
crosses the 1sobaths is given bv the angle 6 . The coordinate system is illustrated in 
Figure 1 on page 10 and includes gravitational. g, and rotational, Q, vectors. 


The equations for the path flow are? 





T D (10) 
dc 

and 
ee (11) 
Che 


The plume is considered to be a well-mixed, broad. thin layer of high density fluid adja- 
cent to the bottom. Pressure gradient forces pull the plume down the slope while 
Coriolis forces bend the flow to the right. If frictional forces were not considered, the 
flow would move horizontally along the isobaths (Killworth. 1977) in geostrophic bal- 
ance. Friction is included, using quadratic drag laws to relate frictional resistance to the 
Square of the mean velocity (Smith, 1975). The plume ts thus allowed to flow down the 
slope, constantly entraining ambient fluid and losing its negative buovancy. Once the 
plume reaches neutral buovancy, the plume interleaves with the ambient waters. The 


equations of the plume’s flow from Killworth (1977), are: 
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LAN = EATY, (12) 
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Figure 1. The streamtube coordinate svstem as given by Smith (1975). 


y? P. = A sin 0 cos $ — fV. (13) 


These equations are, the along-stream derivatives of the mass flux (12), the momentum 
flux (13), the buovancv flux (14) and the cross stream momentum balance (15). In the 
above equations, A is the cross-sectional area of the plume, V 1s the plume’s velocity, A 


is the plumes buovancy, @ 1s the slope angle, f ts the angle at which the plume crosses 


the isobath, N 1s the Brúnt-Váisálá frequency, K, is the drag coefficient. fis the Coriolis 
parameter, and £, is the entrainment coefficient. 


The equations were integrated using a modified Adams method, 


Predictor: 


/ =. = ^ 
a + (16) 


Corrector: 


/ " 
dt ee =T Deo N m). (16a) 


with the initial steps provided with the Runge-Kutta method, 
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This model uses  Killworth's (1977) entrainment and frictional 
parameterizations, which has the entrainment rate and frictional force porportional to 
the product of the area and velocity (velocity squared for friction). Thus, 
PC" E = Es y ‘where w is the width of the plume and b is the depth of the plume. 
Also used in this model is Bo Pedersen’s (1980) entrainment parameterization for small 


Slopes. Melling and Lewis (1982) included this in their version of Smith's model. The 
K 
Ri 





entrainment coefficient is given as £ 20.072 sin 0 sin fj 2 0.072( 
Ri 2 b cos 0| V? 1s the bulk Richardson number of the flow. 
5. The coupled model 


) where 


The model takes the form illustrated in Figure 2 on page 13. Underneath the 
polvna, a plume with a depth, b. and a width. w, 1s salinized due to brine rejection. The 
plume has an initial velocity, V,. The effects of varying the initial velocity will be ex- 


amined in Section 3, however, the standard case value of Y, = .04m/s 1s used. This value 


is close to the values used bv Melling and Lewis (1982) and observed values of currents 
adjacent to polvnas reported bv Schumacher et al. (1985). The value of V, determines 
the areal flux used in calculating the salinity increase in the water column. The initial 
Velocity 15 increased after the non-entraining. non-rotating rudiment plume travels 10 
km straight downslope. Effects of rotation are neglected at this early stage of the plume 
to minimize the exposure to brine rejection under the polvna. Increasing the velocity 
after 10 km ensures more realistic salinity values of the plume, and would, in any case, 
occur in the first steps of the streamtube model. The non-entraining, non-rotating 
rudiment plume continues straight downslope until a depth of 45 meters 1s reached. This 
depth was chosen due to its use bv Melling and Lewis (1982) as the initial depth of the 
streamtube model. 

At 45 meters, the plume descends the slope according to the dvnamics set forth 
and established by Smith (1973), Killworth (1977), and Melling and Lewis (1982). Three 
stratified lavers are defined in this coupled model. Layer 1 lies between 45 and 55 meters 
depth, Laver 2 lies between 55 and 300 meters. and Laver 3 lies below 300 meters. The 
plume is assumed to interleave with ambient waters when the velocity is less than .OL 
m's: however, all sensitivity studies will be carried out using 90 days as the maximum 


number of days the plume can descend while retaining its integrity. 


B. MODEL PARAMETERS 
l. Environmental variables 

[n order to examine the environment's influence upon possible polyna- related 
production of halocline and deep waters, certain environmental variables were individ- 
uallv varied during model simulations. The sensitivity of the model to a variable's al- 
teration 1s defined bv the dependence of the plume's depth after 90 days upon the 
particular variable, relative to a standard case, The relationship between the depth of 
the plume and a single environmental variable will be reported in Chapter 3. [he vari- 


ables to be investigated are listed in Table 1 on page 14. 
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Figure 2. Schematic representation of the coupled model. As the shelfwater is ex- 
posed to the polyna, a high density plume is formed, resulting in the 
plume's descent down the slope until it 15 neutrallv buovant. 
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Table 1. LIST OF ENVIRONMENTAL INPUTS: The dependency of the plume 
depth after 90 davs on each of the variables below will be examined in 
Section > 


air temperature 









wind velocity 


shel G ater SE 


secondarv slope 
Drünt-Váisálà frequency. laver 1 










Brünt-Väisälä — laver 


2. Environmental conditions for the standard case 

To achieve consistency with pre-existing works, manv of the environmetal con- 
ditions for the standard case were borrowed from Melling and Lewis (1982), Killworth 
and Smith (1984), and Pease (1987). The standard case initial values are listed in 
Taolet? oni paretin 

The most significant departure from an initial value given by previous work, is 
the H, value. In this paper, the frazil ice collection thickness is taken to be 22 me 
Pease (1987) shows a higher collection thickness results in longer polyna opening time 
and larger polvna width for the same wind speed and air temperature. Weeks and 
Acklev (CRREL Monograph S2-1, 1982) give higher values for frazil ice collection 
thicknesses in a wind driven scenario vice .01 to .1 meters in quiescent conditions. Ín 
addition to the frazil ice collection thickness, the salinity of the ice 1s set to 7 ppt (Cox 
and Weeks, 1974) instead of 5 ppt used in Killworth and Smith (1984). 


Table 2. INITIAL VALUES USED FOR THE STANDARD CASE. 
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A. STANDARD CASE 

In order to examine the environmental influence on the possible polkna production 
of haloehne and deep waters, a standard case was defined with values which are tvpical 
of Arctic Ocean observations. The offshore wind velocity requires a considerable 
meteorological event providing winds higher than the monthly mean winds. However. 
the value of 15 m’s for offshore wind velocity is within observed wind speeds OBESSE 
storms reported in previous papers (Schumacher, et al. (1983) and Aagaard mem 
(1983)). The standard case plume obtains a salinity of 34.98 ppt after transversing the 
polvna. which corresponds to a density of 1032.93 gm/n? and a negative buoyancy of 
0643 ın/s? . It is therefore not suprising when the plume penetrates the same Arctic 
pyenocline used in Melling and Lewis (1982) to a depth of 436 meters. Figure 3 on page 
17 shows the velocity of the plume as it descends the slope. Due to the plume’s high 
negative buovancy, the plume accelerates slightly on the initial slope then rapidly in- 
creases in speed after the shelf break at 33 meters. The maximum velocity 1s 0.30 m/s 
reached at a depth of 60 meters. The mass flux, AV, increased bv over 1600% empha- 
sizing the dependence of the entrainment rate on the velocity of the plume. At the end 
of 90 days. the standard case plume has a velocity of 0.017 m s and a negative buovancy 
of 0.000482 m/ sec?. The plume has travelled 253 kilometers along the shoreline and 96 
kilometers down slope. The flow path of the plume can be seen in Figure 4 on page 
18. Overall. the plume in the standard case behaves similarly, although with different 
magnitudes, as the plumes in Melling and Lewis (1982). With this in mind, sensitivity 
studies of the impact of individual environmental parameters on the 90-day plume depth 
can be exanuned. Sensitivity experimental results typically include 10-20 data points with 


increased resolution in areas of abrubt change. 


B. SENSITIVITY EXPERIMENTS 

The polyna's maximum width is dependent on the growth rate of ice, the advection 
rate of ice, and the collection thickness of ice. Air temperature and the wind speed af- 
fects the growth rate of ice, while the advection rate 1s proportional to the wind speed. 
Lowering the air temperature results in a monotonic increase in the plume’s depth at 90 
davs as shown in Figure 5 on page 20. This is expected due to the increased production 


of frazıl ice resulting in higher salinization of the water column. Although Pease (1987) 
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Figure 3. The standard case plume”s velocity is shown with depth: Note the rapid 


increase in velocitv as the plume crosses the shelfbreak at the 55 meter 


isobath to a steeper slope. 


has shown the polvna size decreases with lowering air temperature at a given wind speed, 
the plume’s salinity (hence, density) increases most rapidly in the first 10 kilometers of 
downslope travel. Here the plume will be subjected to significantly higher salinization 
due to higher ice growth rate and lower areal flux rate. Although the polvna size is 
smaller with lowering air temperatures, the majority of the salinization has already oc- 


curred in the first 10 kilometers. Indeed. the polyna shows a decreasing salt contribution 
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Figure 4. The standard case plume flow path in the x-y coordinates: The effects 


of Coriolis can be seen as the plume's angle f decreases during the 


plume's descent. 


BITE TIIG EET RET decreasing air temperature. This cam be seen in 
Table 3 on page 19. This decreasing salt contribution to the water column is due to the 


polvna’s smaller width. 


Table 3. SALINITY CONTRIBUTION OF POLYNA WITH DECREASING AIR 
TEMPERATURE: The table shows the inittal salinity of the plume. S, 
EE ae E Kilemeters of downslope travel. the final salinity of the 
piume. S. after travelling the entire width of the plume. x and the de- 
creasing salt contribution of the polvna during the distance A, — lûk with 
decreasing air temperature. 7.» 
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On the other hand, if air temperature is held constant, and offshore wind speed is 
miercased, the polyna grows in size due to increased ice advection rate. The initial con- 
tribution of salt during the first 10 kilometers of the polyna does not varv however, the 
increased polyna size exposes the water column to further salinization. This results in 
a denser, deeper penetrating plume with increasing wind speed, (See Figure 6 on page 
21). The dependence of plume depth upon the wind velocity and air temperature illus- 
trate the possible atmospheric influences upon the Arctic Ocean's salinity characteristics. 
The frequency and severity of Arctic storms which would be conducive to high density 
plume formation is a consideration in the calculation of production rates, but cannot 
be discussed quatitativelv using this simplified model. 

The effect of varving the initial salinitv of the shelfwater on plume depth is shown 
in Figure 7 on page 22. The results of this experiment show deeper penetration with 
higher initial salinitv. This is not suprising due to the significant role of salinity upon 
water density near the freezing point. The plume depth after 90 days is greater than the 
depths reported by Melling and Lewis (1982), even with relatively low salinity values. 

Initial shelfwater velocity defines two parameters within the model, the areal flux 


used in computing the new salinity of the water column and the initial velocity of the 
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Figure 5. Effect of decreasing air temperature on the plume’s 90-day depth: With 


ge Ñ 


decreasing air temperature, the increased ice growth rate results in higher 


initial salinization of the water column. 


plume. The higher the initial velocity, the less saline the plume. Flowever, when using 
the plume depth at 90 days as a standard criteria to examine the results, the higher ve- 
locity plumes reach deeper depths at 90 davs than the denser, but lower velocity plumes. 
If a speed criteria were used, then there would be a negative correlation betweengebs 
plume depth and initial shelfiwater velocity. Instead, Figure 8 on page 23 shows a pos- 


itive correlation between plume depth and initial shelfwater velocity. 
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Figure 6. Effect of varving offshore wind speed on the 90-day plume depth: The 
polvna's increased width with increasing wind speed results in a higher 
density plume duerto longer exposure of the water column to the 


polvna's brine rejection. This results in a deeper penetrating plume. 


Similarly, the parameters of shelf slope and secondary slope appear to have contra- 
dictory results. Figure 9 on page 24 shows as the initial shelf slope is lessened, the 
plume penetrates to deeper depths after 90 davs of travel. In this experiment, the shelf 
width was fixed at 20 kilometers after the 45 meter isobath. This deeper penetration 


appears to be due to a lower entrainment parameter, which is a function of the shelf 
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Figure 7. The effect of initial shelfwater salinity on the 90-day plume depth: With 
increasing shelftrater salinitv the plume attains a higher density due to 


increased salt rejected bv ice growth added to the shelfwater. 


slope. In contrast. Figure 10 on page 25 shows deeper penetration with an increased 


secondary slope. Although the entrainment parameter sull increases with slope. the 
90-dav plume depth 1s limited by the slope itself. Since 90 davs was considered a rea- 
sonable time scale for the plume to maintain its integrity, deeper plume penetrations 


would most likely occur with steeper slopes after the shelf break, ie. in submarine 
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Figure 8. Effect of varving initial shelfwater velocitv on 90-dav plume depth: An 


increased initial Velocity results in a deeper penctraung plume at 90 days. 


canvons. Varving the shelf width after the 45 meter isobath has little effect on the 
90-day plume depth. This is shown in Figure 11 on page 26. 

Examining the stratification of the Arctic Ocean is accomplished by varving the 
Briint-Vdisala frequency of Lavers | and 2. The deep Layer 3 remains constant at a 
value used in Killworth (1977) for deep waters. Laver 1 is the initial pycnocline and re- 
presents the transition from the fresh surface laver to the saltier halocline. Figure 12 


on page 27 shows the 90 day plume depth is little changed even with a substantially 
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Figure 9. This figure shows the effect of the initial shelf slope on the 90-day plume 
depth: Due to the entrainment parameter being a function of slope, a 


steeper slope results in higher entrainment of ambient water. thus more 


rapid loss in negative buovancy. 


weakened initial stratification with M? 2 2.0x10-*sec?. The short length cof tument 


space which the plume is exposed to this layer is probably the cause of its diminished 


impact on the 90-day plume depth. 
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Figure 10. The effect of varving the secondary slope on the 90-day plume 
depth: Even with a higher entrainment rate. the plume reaches a 


deeper depth with a steeper slope. due to the physical boundary of the 


ocean bottom. 


On the contrary, Figure 13 on page 28 shows the dramatuc difference in the 90-dav 
plume depth with small changes in the Brünt-Váisálà frequency in this 245 meter thick 


laver. Additionally, if the second laver's thickness 1s decreased, the plume will penetrate 


10.0 





=O 
= 
he 
= 
QS 
U] 
O o 
U o 
=e 
=] 
as 
= 
SS 
Os 
LOAN 
O w 
Cc» 

e 

o 

Cc 

m 

0.0 1.0 2.0 3.0 4.0 5.0 
SHELFWIDTH (M) x10* 


Figure 11. The effect of varving the shelf width (after the 45 m isobath) on the 


90-day plume depth: Varying the shelf width has little effect on the 


90-day plume depth. 


to deeper depths. This is illustrated in Figure I4 on page 29. Note the discontinuity in 
the slope of the curve 1s caused when the plume no longer pentrates the second laver. 


Next. the standard case is compared in detail with a “submarine canyon” simulation. 
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Figure 12. The effect of varying the initial pycnocline strength on the 90-day plume 
depth: The initial pvenocline. between 45 and 55 meter depths. has a 


small effect on the 90-dav plume depth due to the small thickness of 


this larer. 


C. THE "SUBMARINE CANYON" CASE 
The “submarine canton case is a Modification of the standard case using a steeper 


Secondary slope, lower Layer 2 Drünt-Váisáli frequency value, anda decreased Laver 2 
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Figure 13. The effect of varving Layer 2 Brünt-Väisälä frequency on the 90-day 
plume depth: Due to the plume's longer exposure time in this laver, 
changes ın Brünt-Väisälä frequency tend to have a more dramatic eflect 


on the vertical motion of the plume. 


thickness. [n addition. the angle $, is forced to the initial value of 29 degrees to simulate 
the steering of the plume by topography. The initial values used in the "submarine 


canvon” case are shown in Table 5 on page 31. Forcing the angle f towards its inital 
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Figure 14. Decreasing Layer 2 thickness results in a deeper penetrating 


plume: Xs the plume can no longer penetrate Layer 2, due to its 


thickness. the 90-dav plume depth becomes constant. 


value results in the plume travelling down the Y axis further and less down the X axis. 
Of course, the value chosen for forcing fj will alter the X-Y flow path of the plume. The 
Momepati Of the submarine canyon" case is seen in Figure 15 on page 32. The plume, 
descends to a depth of over 1300 meters in 90 days. has moved 207 kilometers along- 
Shore and 160 kilometers off-shore. In Figure 16 on page 33. the “submarine canyon“ 


plume is seen to reach a maximum velocity 0f 0.36 m.s whereas the standard case plume 


S 


has a maximum velocity of 0.30 m s (Figure 3 on page 17). The effect of increased 
slope, the angle f, and velocity on entrainment is seen when comparing the volume flux 
of the “submarine canyon” plume (Figure 17 on page 34) and the standard case plume 


(Figure IS onspasen 


D. RESULTS USING OBSERVED TEMPERATURE AND SALINITY FIELDS 

The values for this experiment are the sanıe as the standard case except the Brünt- 
Vaisala frequency is computed using observed values of temperature and salinity. The 
temperature and salinitv records are from Ostlund (1987) and Aagaard et al. (1981) 
which were taken in the Canadian Basin. The results of the experiment are shown in 
Table 4. The 90-day plume depth and maximum velocities are Very sinular to the results 
from the Laver 2 Briint-Vaisala frequency sensitivity experiment. 

In an attempt to model a current as observed in Aagaard et al. (1985). a plume 
width of 20 kilometers was used together with a steep secondary slope and the 
temperature salinity fields reported in Aagaard et al. (1981). The plume’s velocity with 
increasing depth is shown in Figure 19 on page 36. Aagaard et al. (1985) reported a 25 
kilometer wide, 15 meter deep plume with a mean velocity of 0.45 m/s. The simulation 
vields a velocitv of 0.30 m.s before the shelf break and 0.92 m,s while descending the 


steep secondarv slope. 


Table 4. THIS TABLE COMPARES SOME RESULTS FROM USING ACTUAL 
T-S RECORDS: The AIWEX-2 and AIWEX-S stations are from 
Ostlund(1987) and the third record is taken from Aagaard et al. (1981). 
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Table 5. INITIAL VALUES USED FOR THE “SUBMARINE CANYON” 
CASE: The highlighted values are those parameters which differ from the 
standard case. 
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Figure 15. The “submarine canyon” plume's flow path: The plume's flow path is 


forced to the initial £ value to simulate the topographic steering within 


a submarine canvon. 
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Figure 16. The velocity of a “submarine canyon” plume with increasing depth: The 


plume has the characteristic "jump" in velocity after the shelf break 


where it reaches a maximum velocitv of 0.56 m.s. 
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Fioure 17. The normalized mass flux of the "submarine canyon” plume: The 
g i p 
plume increases in mass bv over 6800% indicating only 1 5%o of the 


Water reaching the deeper depths 1s shelf water. 
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The normalized mass flux of the standard case plume: 


case plume. 
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Figure 19. Velocity shown with depth of a 20 kilometer wide plume: This figure 


shows a wider plume can penetrate deeper than the narrower 1 kilo- 


meter plumes used in the previous experiments. 
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[his paper has attempted to provide a mechanism in which high salinity, low tem- 
perature water is produced and transported off the continental shelves into the Arctic 
halocline or deeper depths. The environmental parameters which affect the production 
and transport of this cold, salty water were examined, resulting in the conclusion that a 
polvna can produce water with a salinity greater than 34.5 ppt near the freezing point 
which can be transported via gravity flow plumes. The results show deeper penetrations 
to depths greater than 700 meters 1s possible depending on the set of environmental pa- 
rameters. [his indicates polvna-produced high density plumes may also be responsible 
for the ventilation of the deep Arctic Basins, specifically the Canadian Basin. 

Ostlund et al. (1987) placed the residence time of Canadian Basin Deep Water be- 
tween 500 and 800 vears. Using the mixing ratio of 2:1 (shelfwater to intermediate wa- 
ter) from Aagaard et al. (1985), shelf water with a salinity of 35.1 ppt requires an 
off-shelf flux of approximately .0063 Sv to ventilate the Canadian Basin in 800 years. 
The results from the 20 kilometer wide plume experiment give an annualized flux of 
NOS) SV LA Other words, an averase of about two such events per year where 
polvna-produced plumes penetrate below the temperature maximum over the past 800 
vears could account for the ventilation of the Canadian Basin. Furthermore, Ostlund 
et al. (1987) states that the deep water probably consists of much less than 10*o 
shelfwater and would have to be a salinity of 36.15 ppt. The “submarine canvon” case 
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shows the plume reaching a depth of over 1300 meters consisting of onlv 1.3% 
shelfwater at this depth. Table 3 on page 19 shows under extremelv cold conditions, a 
Salinity value of 36.23 ppt can be achieved with this model. Therefore, it is reasonable 
to expect that polvna-produced saline water mav significantly contribute to the venti- 
lation of the Canadian Basin. 

Climatological variations which would affect the number of deep penetrating plumes 
are worthy of discussion. Besides the major periods of global ice ages and their 
counterpart warm eras, there have been recorded several shorter cold and warm epochs. 
These shorter periods coincide with weak and strong wind circulation, respectively. In 
the 20th century, both a warm period lasting about 20 vears and a colder period begin- 
ning around 1950 have been recorded. In each of these two eras, there has been a sud- 


Geneoullow Of surface lavers from the Arctic Ocean. Increased ice flow was recorded 
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in 1958 around Iceland during the warm period. while an increased flow of cold. low 
salinity Water reached Iceland from the Arctic in the 19605. Removal of the colder. 
fresher surface layer from the Arctic Basın would reduce the vertical temperature and 
salinity gradients resulting in a less stable vertical structure for the ocean. In effect. an 
environment simular to reducing the model's Layer 2 Brúnt-Viisálá frequency and re- 
ducing the Laver 2 thickness may occur. enhancing deep plume penetrations. In addi- 
tion to possible weakening of the vertical structure, removal of the ice cover leaves large 
areas of open water where heat is more readily removed from the ocean. In these large 
areas of open water, frazil ice grows rapidlv and abundantly to increase the brine rejected 
into the water column. [n the warm period of the early 20th centurv, the extent of the 
Arctic's semi-permanent ice pack was greatly reduced (up to 10%) between 1920 and 
1938. The reduction in the semi-permanent ice pack is attributed to higher ice removal 
rates (increased advection) due to the strong wind circulation of the period. In Ostlund 
et al. (1987). it 15 noted that the deepest waters in the Canadian Basin may be the re- 
mains of the warm climate which occurred around 1000 A.D. vice any cold epochs such 
as the "Little [Ice Age" in the late Tth and cari Lohi entunes 

In summary. periods when Arctic wind circulation is strengthened resulting in open 
water due to rapid ice divergence coupled with a weaker vertical Arctic Ocean sub, 
seems to be the most conducive environmental scenario for deep plume penetrations to 
provide the additional salt source in the Canadian Basin. 

In Killworth and Smith (1984), the heat diffusivity constant, « , required to achieve 
steady state balances requires some peculiar changes in magnitude with increasing depth. 
A value of k=3 x 10-5m?s*! 1s established to achieve a steady balance between the non- 
turbulent shelfwater plume used in Killworth and Smith (1984) and the region where the 
plume interleaves (note: the densitv of this plume is based on the brine rejection due to 
a single vear's ice growth over a continental shelf of 50 meters average depth). To 
properly supply the East Greenland current outflow, xk was required to be 
=3 x 10m? secr!, And finally, to balance the inflow of warm Atlantic water and upward 
advection of heat, « at this level was required to be =1.2 x 10 “wn? sec If the nomi 
mass flux of the turbulent plume model was to be used as a gauge to understand the 
turbulent diffusion, this thesis would present some interesting results. In Melling and 
Lewis (1982), the normalized mass flux was of order =10%, In the standard cas 
normalized mass flux was of order =10!. The “submarine canyon” case had a mass flux 


of order =10?%, Could it be that the three above cases provide a insight into the missing 


phvsics so stated ın Kıllworth and Smith (1984)? Some other missing physics may be the 
double diffusion process as outlined in Carmack and Aagaard (1973) and the breaking 
of internal waves (Perkin and Lewis, 1978). Aagaard et al. (1985) suggests the incorpo- 
ration of turbulent plume dynamıcs into the “filling box” model of Killworth and Smith 
(1984). A mechanism providing a specific low temperature, high salinity water source 
has been developed in this paper. Although the mechanism is episodic. as suggested 
Aagaard et al. (1985), (Killworth and Smith (1984) had earlier used the term 
“spasmodic ), due to the dependence on environmental variables, flux rates of this highlv 
saline water can be estimated from renewal times of the Canadian Basin. The problem 
of incorporating the turbulent plume in the “filling box” model mav be tackled, but is 
bevond the scope of this paper. 

The maintenance of the Arctic halocline by the mechanism presented in this studv 
alone would require many ümes the number of active polvnas producing saline 
shelfwater. Using the 3:2 nuxing ratio proposed by Aagaard et al. (198i. a minimum 
of 2.5 Sv production rate of cold. saline shelfwater was determined necessary to replenish 
the halocline over a period of 10 vears. With the flux rate of the 20 km wide plume (.0035 
Sv), the individual mechanism modelled in this paper would require 200-300 coastal 
polvnas producing 2-3 high salinity, low temperature plumes per vear to maintain the 
Arctic halocline. Although 200-300 such polynas may be thought prohibitive. areas of 
open water which exist during the onset of freezing temperatures in the late fall, as well 
as leads which continue through winter over as much as 10% of the ¡ce pack, could 
provide a rapid source of salt to help overcome the summertime salinitv gradient. In 
times of strong vertical structure, these polyna-produced plumes can transport low 
temperature, high salinitv water into the Arctic halocline. Overall, the polyna seems to 
be an important source of brine which may help maintain the Arctic halocline and ven- 


tilate the deep basins. 
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